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The properties of bosonic Ytterbium photoassociation spectra near the intercombination transition ^So-'^-Pi 
are studied theoretically at ultra low temperatures. We demonstrate how the shapes and intensities of rotational 
components of optical Feshbach resonances are affected by mass tuning of the scattering properties of the 
two colliding ground state atoms. Particular attention is given to the relationship between the magnitude of 
the scattering length and the occurrence of shape resonances in higher partial waves of the van der Waals 
system. We develop a mass scaled model of the excited state potential that represents the experimental data for 
different isotopes. The shape of the rotational photoassociation spectrum for various bosonic Yb isotopes can 
be qualitatively different. 
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PACS numbers: 34.50.Rk, 34.10.+X, 34.20.Cf, 32.80.Pj 



I. INTRODUCTION 

The properties of intercombination transitions in alkaline 
earth atoms and atoms with similar electronic structure has 
become an object of a growing number of experimental and 
theoretical studies. It is mostly caused by a variety of new 
applications in the physics of ultra cold atoms: from laser 
cooling and trapping jll 0, [U to optical frequency standards 
61,151,161,11. 

Great progress in this area has been achieved for Ytterbium 
(Yb), which has 7 stable isotopes with atomic weights 168, 
170, 171, 172, 173, 174, and 176. Quantum degenerate gases 
have been obtained for the bosonic isotopes ^^"^Yb HI], ^^"Yb 
il], and ^'''^Yb ll^l, and the fermionic isotopes ^"^^Yb and 
i73y1j |[Yi1i , Photoassociation spectroscopy has been carried 
out for bosons lfl2l[l3ll as well as fermions ifTill . A two-color 
photoassociation experiment has allowed a precise determina- 
tion of the ground state s-wave scattering lengths for all com- 
binations of Yb isotopes |jT5|. Finally, it was demonstrated 
that the significant change of the scattering properties for col- 
liding ground state atoms can be achieved with optical Fesh- 
bach resonances in these systems lfl6ll . 

This work is devoted to the theoretical study of the photoas- 
sociation spectra near the intercombination transition ^5*0- 
■^Pi of Yb for bosonic isotopes. We take advantage of the 
experimental photoassociation spectra for ^^"'Yb2 and ^^^Yb2 
from Tojo et al. || 1 311 to precisely determine the binding ener- 
gies of excited state molecular energy levels. In addition, we 
report new measurements of the binding energies of excited 
170 and ^^^Yb2. Analysis of the photoassociation spec- 
tra also requires knowledge of the scattering properties in the 
ground electronic state of two colliding atoms. The necessary 
information and experimental data for all ground state isotopic 



combinations can be found in Ref. ifTsll . 

Yb is an excellent example of a system for which the scat- 
tering properties can be easily tuned by the change of the iso- 
topic combination, thus changing the reduced mass of the col- 
liding pair. Such mass tuning of the scattering properties can 
also be very useful in other similar species with several iso- 
topes. Development in the laser trapping and cooling of atoms 
such as Cd flT] and Hg iHKll will hopefully allow photoas- 
sociation investigations of these systems in the future. Cad- 
mium and Hg, like Yb, are good candidates for mass tuning 
of the scattering length because of their numerous isotopes. 
In addition, Hg is seen as a very promising candidate for fu- 
ture optical frequency standards. The clock frequency shift in- 
duced by black body radiation is especially small in Hg lfl9ll . 
compared with other Group II elements 1I20I1 . 



II. PHOTOASSOCIATION RESONANCE 

The two-body loss rate coefficient K{A; I, T) in the pho- 
toassociation process for a thermal cloud of ultracold atoms 
at temperature T needs to be evaluated as an average over all 
possible momenta of two colliding atoms. This loss rate is 
directly dependent on the light intensity / leading to the pho- 
toassociation as well as on the detuning of the light from the 
atomic resonance A and detuning Ag corresponding to the op- 
tical resonance coupling the scattering ground state "5" with 
the excited "e" bound state fl^. 

The averaged loss rate can be written as ll22ll : 



K{A;I,T) 



(/C(A;/,p;,p-))- 



(1) 



where 
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/C(A;/,p„p;) = 1^ ^(2J, + 1) — 



(2) 



r 



is the loss rate II23L 1241 12511 corresponding to particular mo- 
mentum vectors of the relative motion of the two colliding 
atoms pr as well as the motion of their center of mass pc- Con- 
tributions from all possible transitions between excited bound 
and ground scattering states are included in this expression. 
They are taken in the sum with weights dependent on the to- 
tal angular momentum Jg of the two-atom system. In Eq. ^ 
the magnitude of the wave vector corresponding to the relative 
motion is kr = Pr/h, the kinetic energy of relative motion is 
Sr = ?l^fc^/(2/i), and fi is the reduced mass of the colliding 
atoms. The Doppler shift is described by ££) = — fikiiis-PcJ-M 
where the magnitude of the laser light wave vector is fcias = 
uj/c and the mass of the molecule created in the photoasso- 
ciation process is Ai. The shift of the photoassociation res- 
onance £"100,11101 = ^^^fas/(2A^) caused by the photon re- 
coil is also included here. Finally, the light induced shift of a 
given photoassociation resonance can be expressed as a sum 
Ee{I, £r) = J2g ^eg{I, £r) of Contributions Eeg{I, Er) of all 
possible optical transitions between the excited bound state 
"e" and ground scattering states "g". Similarly, the total width 
of the resonance Te{I,er) = Fpe + J2g^eg{I,£r), where 
Feg(/, Er) is the light induced width between the "e" and "g" 
states and Fpe describes all possible other processes leading 
to loss of the excited state. If radiative decay is the dominant 
loss process, then Fpe can be taken as the natural width of the 
excited molecular state. 

The light induced width ||21|24 



Teg{I,er) = 27r Vus{I) (er))| 



(3) 



is linearly dependent on the light intensity / through the op- 
erator Vias (/) describing optical coupling between particular 
excited and ground states. For the case investigated here more 
details about this operator can be found in Refs. i23,i26ll . This 



width also depends on the kinetic energy of the relative motion 
of the two colliding atoms This is because the energy nor- 
malized ground scattering state | ^1'+ (e^ ) ) strongly depends on 
this energy. Finally the magnitude of the light-induced width 
is dependent on the unit normaUzed excited bound state |5'e). 

The light induced shift ll25ll can be calculated from the Fano 
theory liZAi28il 



Eeg{I, £r) 



1 

2^' 



de 



E 



\{^,\V,,s{I)\^g 
Er- Eg 



(4) 



where "P is a principal part integral over all collision ener- 
gies, and the sum occurring in this expression is taken over all 
unity normalized bound state of the ground electronic 
potential. Like in the case of the light induced width, the light 
induced shift is linearly dependent on the laser intensity / and 
of course on Sr- 

The photoassociation process can be viewed as a case of 
the optical Feshbach resonance ll25ll29l[30ll . The optical cou- 
pling of the ground scattering state to the excited bound state 
changes both the amplitude and the phase of the scattering 
wave function. The amplitude can be detected through loss of 
atoms from the trap. The phase change can also be detected 
lfT6l[3lll . This approach leads to conclusion that the scattering 
length describing ultra cold collisions of the two ground state 
atoms can be modified by light due to the coupling with an ex- 
cited bound state. This effect was demonstrated with a Bose- 
Einstein condensate (BEC) of ^^Rb ll32l[33ll . The expression 
for a scattering length modified by light can be written in the 
following form ll25ll : 



J 



a(A,/) = flbg 



Fpe [A - Ae - Eejl, 0) - ^icc.mol] 
[A - Ae - E,{I, 0) - £,cc,mol]2 + [re(/, 0)/2f 



(5) 



r 



where abg is the background scattering length of the system 
in absence of the Ught. The optical length Er) ll30ll34ll 



reg(/,g.) 
2krT 



(6) 



Because of the s-wave threshold law, this quantity approaches 
an energy-independent constant as ^ that varies linearly 



with light intensity /. The optical length characterizes the 
strength of an optical Feshbach resonance, namely the abil- 
ity of light to change the scattering length of the ground state 
atoms; see also Ref. llssll for a review of Feshbach resonances, 
including optical ones. The scattering length can be changed 
on the order of its background magnitude jobgl while mini- 
mizing losses if the optical length is very large compared with 
I flbg I so that the optical control can be achieved at large de- 
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tuning. 



III. DERIVATION OF A SINGLE CHANNEL MODEL 

The particular properties of the photoassociation reso- 
nances are strongly dependent on the properties of the atomic 
interaction of the colliding atoms. In general we describe 
the colliding system using the Hamiltonian operator H = 
T + Ha + Vint + ^ot- In this expression T is the kinetic 
energy operator for relative radial motion. Ha is the atomic 
Hamiltonian operator representing internal atomic degrees of 
freedom, Vint is the interaction operator described by nonrel- 
ativistic molecular Born-Oppenheimer potentials, and V^ot is 
the rotational energy operator. Reference ll22ll gives more de- 
tails for the case investigated here. 

Let us first focus on the interaction in the excited state. In 
this paper we limit our discussion to photoassociation near 
the intercombination transition ^Sq-^Pi. Therefore we will 
discuss the atomic interaction properties only near the disso- 
ciation limit ^Sn +^ Pi . In this case it is convenient to use 



the \jlJM\p) basis. Here j is the total electron angular mo- 
mentum, I is the rotational angular momentum, and J ^ j + 1 
is the total angular momentum. The projections of ,/ on a 
space-fixed z axis is M. Finally p is the total parity. We add 
an index e to all quantities introduced here to indicate that 
they correspond to the excited electronic state. It should be 
noted that Je as well as A/^ are good quantum numbers of the 
Hamiltonian H. In our case with je = 1 and pe = — 1 we 
can solve the Schrodinger equation using only two channels 
for each J^: one with 1^ = — 1 and one with = + 1. 
In this basis the interaction operator V^nt is not diagonal. Its 
matrix elements can be expressed in terms of potentials Vo{r) 
and Vi (r) which correspond to states with 17 = and ~ 1, 
respectively, where is the projection of the total electron 
angular momentum je along the interatomic axis. The other 
components of the Hamiltonian operator, T, Ha, and Vrot, are 
diagonal in our basis. 



The matrix elements of the sum Vint 
in a compact form ll36ll : 



Vrnt can be written 



J 



If-. jf p_ 1 



2XTT^oM + ^V^ir) + B{r)MJe - 1) 



^iMr) + ItItt^iM + B{r){Je + 2)(Je + 1) / ^ + I 



le = Je + 1 



4g±^[Fi(r)-Fo(r)] 



(7) 



r 



where i3(r) ~ h^/{2pr'^). Henceforth we will omit explicit diagonalization of this matrix gives the following eigenvalues: 
indication of the r-dependence of Vq, Vi and B. Analytic 



V. = l\yo^y^Hv.~vJi-^^ + '-^^+2B[MJe+i)+i] 



(8) 
(9) 



r 



These eigenvalues Vq and Vi can be used as effective poten- 
tials for single channel calculations of excited bound states. 
Such an approach allows one to approximate the influence of 
Coriolis coupling on the molecular structure in simple single 
channel calculations. 

The quantum number O becomes a good one in the range 



of r where |Vi — Vo| 3> B, i e., when the anisotropy of the 
interaction between two atoms is much bigger than the rota- 
tion energy. In such a case the eigenvalues and corresponding 
eigenvectors of the matrix in Eq. (jTji can be written in the 
following form: 
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Vo = l/o+B[Je(Je + l)+2] 



Vl - Vl+BJe{Je + l) 



-V^e/(2Je + l) 

v/( Je + l)/(2Je + 1) 
V(Je + l)/(2Je + 1) 

V Je/(2Je + 1) 



(10) 



(11) 



r 



In many cases these potentials provide a successful ap- 
proximation for single channel calculations of bound states 
of real systems. Such an approach is particularly good for 
Yb, because the potential curves at long range are domi- 
nated by the resonant dipole interaction. In the case of Yb 
1^1 ~ ^1 = fCs/^^ 3> B up to very large r. This is a quite 



different case from Sr OTil . where the C3 value is more than 
20 times smaller than for Yb. In the case of Sr the competition 
between the van der Waals and resonance interactions is im- 
portant. In the range of r where the interaction is dominated 
by the resonant dipole term the ungerade potentials for = 
and il = 1 are attractive and repulsive, respectively. Therefore 
only the potential Vo will support a series of bound states hav- 
ing 0+ symmetry near the dissociation limit. The Vi potential 
becomes attractive at small r due to chemical bonding and 
may support a bound state near the dissociation limit. In the 
rest of the paper we will describe excited bound states of 0+ 
symmetry using single channel solutions of the Schrodinger 
equation. 

This simple approach to the description of atomic interac- 
tions allows us to write the appropriate expressions for the 
light induced width and shift: 



same l^- Therefore the light induced width, Eq. (O, is pro- 
portional to the fraction of the total wave function of the ex- 
cited bound state in the channel with 1^ = Ig. This fraction 
can be found from Eqs. ( fTOl i and ( fTTT i. Moreover, the light in- 
duce width is proportional to the square of the proper Clebsch- 
Gordan coefficient \ {JgKMgq\Jf^Me)\, where K = 1 for the 
dipole transition and q = —1, 0, or +1 describes the polar- 
ization of light. This way one can introduce the rotational line 
strength factor 



flfiMg^q) 



2J. 



1 



2J + 1 

' \{jglOO\JeO)f\{jgKMgq\j,AQf , 

(15) 

which is dependent on the light polarization q and the pro- 
jection Mg of the total angular momentum Jg in the ground 
electronic state, see Ref . 113911 . 

The rotational line strength factor given by Eq. (fl4] i can be 
obtained as an average 



JTOt 

J eg ^ 



1 



2Ja 



Ma 



Tg{Mg,q) 



(16) 



and 



over all possible orientations of the total angular momentum 
(12) J^; compare Ref s. Issl [39I1 . This quantity is independent of 
the light polarization q. To derive Eq. ( fT4l ) one can use the 
following relation 



Q r\3 

E,g[I,er)^TA^^f, 



J e 



eg J eg 



Her). 



(13) 



In Eqs. (fTZt and ( fTSl l Ta is the natural decay width of the 
atomic transition and \a the wavelength of the atomic tran- 
sition. The dimensionless rotational line strength factor for a 
transition from the ground scattering state to the excited 0+ 
state in a bosonic isotope has the simple approximate form 
obtained by Machholm ef fl/. i38<l : 



(•rot 

leg — 



1 2Je + l Je 

3 2,/g + l 2,/e + l 

1 2Je + l Je+l 

3 2J„ + 1 2Je + l 



1 J3 + I 

3 2J„ + 1 



3 2J„ + 1 



for Je = Ja 



for Je = Ja 



1 , 



(14) 

for even Jg. This expression can be obtained using the ap- 



proach described in Ref. 112211 . To extract the rotational line 
strength factor /'g*" one can start from Eq. (O and use Eq. 
(B3) from Ref. 1I22I1 . In our case, this allows us to show that 
the laser radiation couples the ground scattering channel hav- 



ing 



Jg only with the excited state channel having the 



Y,\{JaKMgq\JMf = 



2X 



1 



(17) 



which is fulfilled for Jg = Jg — 1 and Jg + 1, A' = 1 and 
Me = Mg + q. Equation ( fT4l i can be also expressed in terms 
of a Wigner 3-j symbol, see Ref. ll39ll : 



l-TC 

J en 



|(JglOO|JeO)f 



2Je + l 



Jg 1 Je 

000 



(18) 



This work is limited only to the case of weak interaction with 
light. Therefore we can ignore the dependence of the light 
induced width and shift on the projection Mg and use the pro- 
jection independent Eq. (fl4b . 

The total angular momentum Jg in the ground scattering 
state is the same as the rotational angular momentum Ig. This 
is a consequence of the fact that the atomic interaction in the 
ground electronic state has ^ S+ symmetry and the total elec- 
tronic angular momentum jg = 0. In homonuclear pairs of 
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bosonic atoms only even Jg are allowed. Finally, the Franck- 
Condon factors 12511 : 



dr 4>e{r)fg{r;£r) 



(19) 



and 



dr (j)e[r)gg[r;er) / dr' <j)e[r')fg{r' -,6,-) 
Jo 

(20) 

can be expressed in terms of the regular fg{r; e,.) and irregu- 
lar gg{r; Sr) solutions of the the Schrodinger equation for the 
ground scattering state and the wave function 4)e{r) for the 
excited bound state. 

The expressions for the light induced width and shift can be 
further simplified by using the reflection approximation ll24l 
The Franck-Condon factors can then be written in the 
following form; 



Fc-r 

eg 



dE, 1 
dn Dc 



\fg{rc;er)Y 



and 



j-FC-E/ \ ^ dEf. 1 



2 dn D 



(21) 



(22) 



c 



where Dc = 9Ve/57'|^^ - dVg/dr\^,^ the fraction ^ is 
the mean vibrational spacing, Vp and Vg are the effective po- 
tentials in the excited and ground electronic states. Here rc is 
the Condon point. When the reflection approximation is ap- 
plicable, the Condon point is approximately the classical outer 
turning point for the excited bound state. 

The reflection approximation is very good for a number of 
Yb bound states close to the ^Sq+'^Pi dissociation limit. This 
is in contrast with other species such as Ca or Sr, for which 
the reflection approximation needs to be modified 1411!]. Since 
Eqs. (fT2] | and (l2Tl i show that the light induced width is propor- 
tional to \fg{rc)\'^, photoassociation spectroscopy provides 
an excellent tool to investigate the properties of the ground 
state scattering wave function 14211 . This will be explored be- 
low for photoassociation of Yb near the intercombination line. 



IV. EXPERIMENTAL DATA 

Experimental information is crucial for determining the pa- 
rameters that describe the long range interaction near the ^Sq- 
•^Pi dissociation limit. We used the experimental data ob- 
tained by Tojo et al. ifisll to find the parameters needed to 
describe Vb(r) for the 0+ state. As described in Ref. ifisll . 
Yb atoms were decelerated by a Zeeman-slowing laser for the 
^So-^Pi transition, and were collected in a magneto-optical 
trap (MOT) with a laser for the ^Sq-^Pi transition. After the 
compression of the MOT, the atoms were transferred into a 
crossed optical trap with laser beams at 532 nm. The atoms 
were evaporatively cooled by decreasing the trap potential 
depth. Typically, MOT duration time is 10 s and evaporative 



cooling time is 6 s. A laser beam for photoassociation was ap- 
plied to the trapped atoms, and the number of remaining atoms 
was measured through an absorption image with the ^So-^Pi 
transition after the release from the optical trap. The applied 
photoassociation light intensity varied between 6.5 /iW/cm^ 
to 90 mW/cm^, depending on which excited level was be- 
ing probed. We obtained the atom-loss spectra by scanning 
the photoassociation laser frequency. These spectra allowed 
the determination of the bound states energies of a series of 
Je = 1 and Je = 3 0+ levels for each of the two isotopic 
species. Table U lists the measured energies along with their 
error bars. The data for the ^^^Yb2 Jg = 1 states were taken 
at about 4 /iK, and the other data in Table U were taken at tem- 
peratures in the range from 5 to 27 /xK. The shift with temper- 
ature and light intensity of the photoassociation features was 
taken into account in the data analysis. These shifts are mostly 
responsible for the magnitude of the error bars listed in Table 

m 

The same technique was used to determine bound states en- 
ergies of the excited homonuclear molecules made from two 
other bosonic isotopes, ^^°Yb and Yb. Long MOT time of 
about 60 s was needed for ^^"Yb because of its small natural 
abundance and the evaporative cooling had to be done in 
a short time of about 1.5 s for ^^^Yb because of three-body 
recombination atom loss due to its large negative scattering 
length liT6il . Table HH lists the values of measured binding en- 
ergies in the excited 0+ Jg = 1 levels of ^^*'Yb2 and ^^^Yb2 
molecules. These data were taken at about 2 /iK. The uncer- 
tainties shown in Table II are mostly due to the light-induced 
shift. 



V. MODEL POTENTIALS 

The location of a photoassociation resonance as the laser 
frequency is scanned is directly related to the binding energy 
of the excited molecule. Ab initio potential curves are still 
insufficiently accurate to describe the position of the most 
weakly bound states. Therefore we introduce an analytic 
model potential valid at large r and fit its parameters to match 
measured bound state energies near the dissociation limit for 
different isotopes. The basic concept is similar to that used 
very successfully for bound state and scattering properties of 
the ground state |[l5il43ll . The potential is chosen to have the 
correct long range form and an arbitrary short range form that 
permits us to represent the correct absolute phase due to the 
unknown short range interactions. This allows us to success- 
fully mass scale the excited state binding energies for different 
isotopic combinations. Mass scahng is more fully described 
in the next Section. 

The interaction potential of two atoms in the 0+ excited 
electronic state for large interatomic separations can be well 
approximated by the following expression (see Ref. IstIi ): 



c. 



(e) 



(23) 
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TABLE I: Comparison of measured binding energies in excited state for ^^''Yb2 and ^'®Yb2 with binding energies calculated from the 
optimal model obtained by least squares fit of experimental data for ^^''Yb2 and ^^''Yb2 isotopes. See details in the text. All quantities are 
given in MHz. 
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-2.4 


-619.3(3.0) 


-616.6 -2.7 


-544.1(5.0) 


-534.1 


-10 


-906.3(3.0) 


-903.9 -2.4 


-798.1(5.0) 


-787.9 


-10.2 




-1296.9 


-1145.1(5.0) 


-1136.6 


-8.5 




-1826.5 




-1608.1 





where the resonant dipole-dipole interaction coefficient 

and (T^^^ is a free parameter that allows us to adjust the phase 
associated with the short range potential. The wavelength 
of light in vacuum corresponding to the transition between 
atomic states ^Sq and ^Pi of Yb is A = 555.802 nm, and 



r is the atomic lifetime. We use the least-squares method 
to optimize the values of cl^\ and cr^'^) while 

matching the calculated binding energies to the experimental 
values measured by Tojo et al. ifisll for ^^"'Yb and ^^^Yb. 
We obtain C^""^ = 0.1949(11) Etal which corresponds 
to T = 869.6(4.5) ns, cf^ = 2.41(0.22) x 10^ Etal, 

C^"^ = 2.3(1.6) X IQS Ehal, and a^""^ = 8.5(1.0) ao, where 
ao w 0.05292 nm and ~ 4.360 x IQ-^^ J. The er- 
rors quoted give the one standard deviation statistical fitting 
error, and do not reflect any systematic errors in the model. 
Adding the Cg' ' coefficient to the model was needed to im- 
prove the quality of the fit for Je = 3 levels. Introducing 
Cg'^'' allowed us to determine the sensitivity of Cg*^' to vari- 
ation of the shorter range of the potential. This sensitivity 
also contributes to the standard deviation of Cg*^' . For com- 
pleteness we give the fitting parameters to enough significant 
digits to reproduce the calculated values to 10 kHz, as shown 
in Table m C^''^ = 0.19488626 Eha^ (r = 869.64762 ns), 
C^"^ = 2405.3647 Ena^, C^*"^ = 229451.31 Ehaf,, and 
cr(^) = 8.4897163 ao. 

Table U compares the binding energies predicted by the 
model to the experimental data for ^^"'Yb2 and ^^^Yb2. The 
model predicts the experimental values to within about one 
MHz on the average, consistent with the experimental error 
bars. 

Having the data for two different isotopes allows us to con- 
struct an appropriately mass-scaled model for the excited state 
levels, similar to what is possible for the ground state ifTsll . 
While cautioning that the short range physics is complicated 
by the interaction of other molecular states with the 0+ state 
and may not be fully represented by a single potential, our 
mass scaled model determines that the most deeply bound ob- 



7 



TABLE III: Corrections to binding energies of the 0„ bound states 
with Je = 1 in the ^^*Yb2 molecule caused by the Coriolis coupling 
and retardation effect; see the text for details. All calculated values 
are given in MHz. 



Binding 




Corrections 




Energy 


Coriolis 


Retardation 1 


Retardation 2 


-4.2 


-0.004 


-0.5 


-0.5 


-9.7 


-0.005 


-0.6 


-0.6 


-20.1 


-0.006 


-0.8 


-0.7 


-38.5 


-0.008 


-1.0 


-0.9 


-69.1 


-0.009 


-1.2 


-1.0 


-117.8 


-0.011 


-1.4 


-1.1 


-192.3 


-0.013 


-1.7 


-1.2 


-302.5 


-0.015 


-1.9 


-1.2 


-461.3 


-0.017 


-2.2 


-1.2 


-684.9 


-0.020 


-2.5 


-1.1 


-993.7 


-0.022 


-2.8 


-0.9 


-1413.0 


-0.025 


-3.1 


-0.5 


-1973.9 


-0.027 


-3.5 


0.0 



served levels at -1973.5 MHz for "'^Yba and -1240 MHz 
for ^^^Yb2 respectively correspond to the v = 106 and 108 
vibrational levels of the model potential. 

An excellent test of mass scaling is to test it using other iso- 
topic combinations. Table HIl shows a comparison of our new 
measured binding energies for ^^"Yb2 and ^^^Yb2 with those 
predicted by our mass scaled single-potential model. There is 
reasonable agreement of about 3 MHz between measured and 
calculated levels, except for the three most deeply bound lev- 
els observed for ^^^Yb2. While 3 MHz is on the order of the 
experimental uncertainty, the accuracy of mass scaling may be 
more limited for the excited states than for the ground states, 
where it was found to be good to approximately 0. 1 MHz lITsIl 
for binding energies up to 325 MHz. Another source of error 
in our single channel excited state model could be the neglect 
of interactions with short range eigenstates of other symme- 
tries. For example, it may be that bound states of 1„ symmetry 
near the threshold perturb the 0+ bound states that are near in 
energy to them. This could be the reason for relatively large 
deviations for the three most deeply bound states of ^^^Yb. 

Our model potential describing the interaction in the elec- 
tronic excited state Oj has a very different shape from those 
obtained in an ab initio calculation by Wang and Dolg liill . 
Clearly such a model can not be treated as a good representa- 
tion of the short range interaction. Nevertheless, its applica- 
bility over a range of isotope masses gives us confidence that 
the number of bound states determined from our model is cor- 
rect to within one or two bound states. The vibronic quantum 
number for the ^^^^Yb bound state at -4.4 MHz is 1 18 with the 
ground state labeled as zero. Also the long range interaction 
should be well described by the model used here. 

We have tested the possible influence of the Coriolis cou- 
pling on the results of our calculations. We have compared 
our calculation used for the fits in which effective potential 



TABLE IV: Comparison of experimental ground state binding ener- 
gies IT^ with the present model. All values are given in MHz. The 
last column is the difference between the experimental values and the 
present model. 



Isotope 


V 


J9 


Experiment 


Theory 


Theory 


Diff. 










[151 Present work 




176Yb 


1 





-70.404 


-70.405 


-70.378 


-0.026 




1 


2 


-37.142 


-37.118 


-37.093 


-0.049 


174Yb 


1 





-10.612 


-10.642 


-10.629 


0.018 




2 





-325.607 


-325.607 


-325.602 


-0.005 




1 


2 


-268.575 


-268.576 


-268.570 


-0.005 


173 Yb 


1 





-1.539 


-1.613 


-1.609 


0.070 


172 Yb 


1 





-123.269 


-123.349 


-123.321 


0.052 




1 


2 


-81.786 


-81.879 


-81.851 


0.065 


171 Yb 


1 





-64.418 


-64.548 


-64.522 


0.104 




1 


2 


-31.302 


-31.392 


-31.367 


0.065 


170 Yb 


1 





-27.661 


-27.755 


-27.735 


0.074 




1 


2 


-3.651 


-3.683 


-3.667 


0.016 



was given by Eq. dTol i with those obtained using the effective 
potential in the form of Eq. ([8]), and find only a small differ- 
ence in binding energies for 0+ bound states with Jg = 1 in 
the ^^^Yb2 molecule. Table Hill shows the magnitude of the 
Coriolis corrections in the column labeled "Coriolis". These 
corrections calculated from the single channel model, Eq. ([8]), 
are below 30 kHz and much less than the experimental error 
bars. The relative importance of this correction increases as 
the bound state energy approaches the threshold. The Coriolis 
coupling mostly play a mar gina l role in the Yb2 excited state 
system. By contrast, Ref. OTT showed it needs to be taken 
into account to correctly calculate two most weakly bound 0+ 
states with Jp = 1 in the ^*Sr2 molecule. 

Similar tests were carried out to check the possible in- 
fluence of retardation effects ll45i l46ll . To do this we have 
replaced in the effective potential, Eqs. ( fTOl i and ( |23] ), the 
standard term describing the resonance interaction —C'^'^/r^ 
by the term in which retardation is taken into account: 

-j_C^/r3)[cos(r/A) + (r/A) siii(r/A)], where A = A/(27r) 
OSl [4611 . The corrections to the binding energies caused by 
this change are listed in the column "Retardation 1" of Table 
Hill These corrections are on the order of a few MHz. In con- 
trast, the same corrections for ^^Sr2 molecule IjtIi are more 
than one order of magnitude smaller. To check whether such 
shifts might be detectable in fitting the data, we have modified 
cr^'^-' so as to change the quantum defect and fit the binding 
energy of the most bound level. The column labeled "Re- 
tardation 2" show the differences from the values calculated 
without retardation. Since these differences are on the order 
of 1 MHz or less, comparable to the experimental uncertainty, 
the present data are not sufficiently accurate to come to any 
definitive conclusions about the observability of retardation 
corrections. 

Reference 03l fits ground state binding energy data for two 
isotopes using a similar form for the ground state potential 
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without the dipolar term. 



C, 



ia) 



ria) 



a 



ia) 



(25) 



This model with mass scaling predicted the observed binding 
energies for 4 other isotopic molecules with an error on the 
order of 0.1 MHz or less. Here we have refitted the ground 
state potential by simultaneously fitting the data for all iso- 
topes in Ref. lUSll . The fit has a slightly better than the pre- 
vious one, but does not represent a significant improvement. 
For the sake of completeness, we give the model parameters 
for the global fit with enough significant digits to reproduce 
the calculated values in Table|IV] C^/^ = 1930.2481 Enal, 



(a) 



194683.32 Eha^ and cr^f) = 9.0240156 



The C^^^ and C, 



(a) 



obtained here as well as in Ref. lo 



agree very nicely with those calculated by Zhang and Dal- 
garno ll47ll . The reported values in Ref. IaT^ are Cg^' = 
2070 EhQ^ and C^^^ = 2.023 x lO^Eha^, respectively, with 
an estimated uncertainty of 10%. It should be emphasized 
that the shape of our potential can significantly differ form the 
real one. Nevertheless our model potential should correctly 
represent the number of vibronic bound states in the ground 
electronic state of Yb2 molecule. Therefore even if our model 
is significantly different from an ab initio potential like that 
reported by Buchachenko et al. ll48ll it should give about the 
same number of vibronic bound states. 

It should be noted that similar analysis to that for Yb ifTsIl 
was carried out recently for the scattering properties in the 
ground electronic state of various isotopes of Sr by Martinez 
et al, using a realistic potential ji^. The scattering properties 
from this work are in very good agreement with results ob- 
tained using a potential derived from Fourier transform spec 
troscopy by Stein et al. 



VI. GROUND STATE SCATTERING WAVE FUNCTION 

As discussed in Sectionlllll the strength and shape of a pho- 
toassociation line is to a large extent determined by the ground 
state scattering wave function at the Condon point for the tran- 
sition. Since the Condon points for photoassociation to near- 
threshold excited states tend to be at quite large internuclear 
separation r, the scattering wave function needs to be known 
only at relatively large r. Consequently, in order to describe 
photoassociation in the ultra low scattering energy regime, the 
detailed knowledge of the atomic interaction at short range 
can be compressed to a very few parameters. 

The most important quantity describing scattering during 
an ultracold collision is the scattering length a. If the long 
range potential has the van der Waals form Vg{r) = —Cq/t^, 
the scattering length is very well approximated by a simple 
analytical relation given by Gribakin and Flambaum ifsill 



a 1 - tan $ 



see also Refs. |5 



(26) 



Here the mean scattering length 



a = 2 '^/^ r(5/4) (2/^C'6/?i^) " is a characteristic length asso- 
ciated with the van der Waals potential, where F is the gamma- 
function. This length also defines a characteristic energy for 
the van der Waals potential, e = ft^ / (2/ia^). The phase $ is 
defined by 



$ = 



^-Vg{r)dr, 



(27) 



where tq is the inner classical turning point of V{r) at zero 
energy. The number of bound states N in the potential is lf52ll 



N = 



1, 



(28) 



where the brackets [. . .] mean the integer part. 

The scattering length varies periodically with phase hav- 
ing a singularity when $/7r = A^ — 3/8. This variation can 
be observed experimentally for different isotopes of the same 
species. If we assume that the interaction potential is the same 
for all isotopes so that mass scaling applies, can be var- 
ied by changing the reduced mass /i. The relation between 
scattering length, the energies of near threshold bound states 
and the reduced masses of the colliding atoms was carefully 
studied by Kitagawa et al. ifTsIl for ground state interactions 
of various Yb isotopes. It is useful to define a reduced mass 
difference A/i needed to change ^{^jl) by tt, that is, to change 
the number of bound states in the potential by one. It is ap- 
proximately 



(29) 



Since N = 72 for ^^^Yb, we see that a mass difference of 5 
atomic mass units is sufficient to change the number of Yb2 
bound states by one. Alternatively, varying fi continuously by 
5 atomic mass units will cause the scattering length to vary 
across a singularity over its full range from — oo to +oo, as 
found by Ref. ifisll . There are actually 7 stable isotopes of Yb, 
and thus 28 different discrete physical values of 2/_i that are 
available in the laboratory using different isotopic combina- 
tions. 

The long-range ground state scattering wave function at 
very low collision energy depends on three basic parameters, 
namely the separation between colliding atoms r, the relative 
kinetic energy of collision Er, and the quantum defect associ- 
ated with the phase $ or, equivalently, the scattering length a. 
The reduced mass /^t is a parameter that allows control of the 
quantum defect. Thus, the long-range scattering wave func- 
tion for the van der Waals system has an universal charac- 
ter that can be expressed by a system-independent function 
of the dimensionless variables r/a, Sr/s and n/Afi. All of 
the needed information about the ground state scattering and 
bound state wave functions can be calculated from a knowl- 
edge of the scaling parameters r/a, Sr/eand {fi— fio) / Afi and 
the value of reduced mass /io for which the scattering length 
is singular in a given system. 

The analytical solutions of the Schrodinger equation are 
known for several class of potentials, thus allowing the deriva- 
tions of compact expressions for the scattering length ifsill . 
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FIG. 1 : (Color online) The lower and upper panels show Vg + Kot 
and Vb + Kot , the respective ground state and excited state potentials. 
The upper curve shows the outer turning points of each calculated 
Ot, Je = 1 bound level for ^^''Yb. The dotted lines indicate the 
experimentally observed range of levels. The arrow indicates the 
mean scattering length d = 75.18 ao. The lower panel shows the 
centrifugal barriers for the long range ground state potential for Jg = 
Ig =0 (s-wave), 2 (d-wave), 4 (g-wave), and 6 (i-wave). An energy 
{Vg + Vrot) /kB = 1 mK is equivalent to {Vg + Vrot)/h = 21 MHz. 



Van der Waals systems such as those discussed here can be 
nicely described and very well understood using analytical 
theory. Gao ifssi l56ll used the framework of quantum de- 
fect theory for a van der Waals system to work out a num- 
ber of practical formulas and results. For example, when the 
s-wave scattering length is singular, there is a bound state at 
E = not only for the s-wave but also for Ig = 4, 8, 12, .... 
On the other hand if the scattering length is equal to a there 
will be zero energy bound states for Ig = 2, 6, 10, .... The 
practical consequence for real potentials is that when the 
s-wave scattering length is near such special values, these 
other threshold bound states show up as shape resonances for 
other partial waves. A shape resonance is a quasibound state 
with enhanced amplitude trapped behind a centrifugal barrier 
that can lead to enhanced photoassociation. While the wave 
functions we show are calculated numerically by solving the 
Schrodinger equation for our potential, the analytic results are 
very helpful for interpreting their features, and for evaluating 
approximations like the reflection formula. 

As we already emphasized before, a photoassociation ex- 
periment can be used to map the square of the scattering 
wave function in the ground electronic state at the Condon 
points corresponding to the different excited bound levels. 
Figure [T] shows on the lower panel the effective potential, 
Vg (r) + B{r) Jg( Jg + l), in the ground electronic state, where 
Jg is the rotational quantum number for the partial wave 
£ = Jg. The Fig. [T] shows the centrifugal barrier for a few 
partial waves with low Jg. The upper panel of this figure 
shows the excited state 0+ potential together with bound states 
and their outer turning points, which are essentially the same 
as their Condon points. This panel shows how the rc value 



that is probed changes with the binding energy of the excited 
bound level. Probing experimental levels with binding ener- 
gies Eb/h between 2000 MHz and 3 MHz allow probing the 
scattering wave function between around 80 ag and 800 ao. 

In order to illustrate some properties of the scattering wave 
function we will begin by setting the collision energy to a con- 
stant value corresponding to Er/ks = 100 /^iK. As we will 
show in the next section, shape resonances could be clearly 
observed at this energy, which is below the top of the g-wave 
centrifugal barrier. Figure [2] shows the dependence of the 
square of the wave function magnitude on /i and r. Chang- 
ing 1^1 corresponds to the change of the quantum defect of the 
colliding system, which changes a and the bound state spec- 
trum. Calculations have been done for a few lowest partial 
waves s, d, g, and i. 

In Fig. |2]the brighter areas show regions of higher ampli- 
tude and the dark lines show the nodes of the wave function. 
The effect of the exclusion of the wave function from short 
range by the centrifugal barrier is especially evident for the 
g-wave and the i-wave. For this Yb system, the singularities 
in the s-wave scattering length occur for 2/i = 167.3, 172.1, 
and 176.9 to„ ifTsll . Near these values, an s-wave node moves 
in to smaller distances on the order of a = 75.18 ao as 2/i in- 
creases and a new bound state appears in the potential. There 
is also a g-wave shape resonance with a dramatic enhance- 
ment of short-range amplitude near these 2/i values where a 
is singular, as expected from the analytic van der Waals quan- 
tum defect theory. Similarly, there are d- and i-wave shape 
resonances near the values of 2/i = 169.7 and 174.5 rriu where 
a = d for the Yb system. Near these shape resonances, a node 
moves to shorter distances inside the centrifugal barrier as 2/j 
increases and there is an extra bound state in the potential. 

If the collision energy £,. is lowered, the wave functions 
show similar patterns, except that the short-range wave func- 
tion is much more attenuated due to lower penetration through 
the centrifugal barrier Clearly, with 28 different physical val- 
ues for 2/i available, we can expect significant isotopic vari- 
ation in the photoassociation spectra, depending on the range 
of Condon points sampled and the temperature of the sam- 
ple. Spectral lines with Condon points outside of the barrier 
will have less isotopic sensitivity, whereas lines with Condon 
points in or inside the region of the barrier will show much 
more sensitivity to the isotopic combination. 



VII. ISOTOPIC VARIATION OF PHOTOASSOCIATION 
SPECTRA 

In this section we will show how the scattering properties of 
the different Yb isotopes affect their photoassociation spectra. 
As an example we performed calculations of the light induced 
trap-loss coefficient at various gas temperatures in range from 
10 yuK, to 1 mK. Results for T =- 100 /iK are particularly in- 
teresting, since the line shapes are relatively sharp and ground 
scattering resonances can have significant influence on the 
spectrum. To emphasize the qualitative differences between 
spectra of the different isotopes we first focus our attention 
on relatively deep bound states having binding energy around 
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FIG. 2: (Color online) The squared magnitude of the numerically calculated scattering wave function \ fg (r, Er) P of two colliding ground state 
Yb atoms is shown as a function of intermolecular separation r and reduced mass of colliding system ^, where m„ ~ 1.6605 x 10~^^ kg . 



The calculations were carried out for energy of collision i 



100 fiK and for Jg = Ig — (s-wave), 2 (d-wave), 4 (g-wave), and 6 (i-wave). 



The dark lines show the nodal lines of the wave functions, whereas the bright areas show the largest amplitudes. 



2000 MHz. As can be seen in Fig. [T]this corresponds to the 
Condon point placed near or inside the centrifugal barriers. 

We start our discussion from ^^^Yb which has a large neg- 
ative scattering length near an s-wave singularity ifTsll . Con- 
sequently, we expect a g-wave shape resonance at low scat- 
tering energy. The upper panel of Fig. [3] shows the energy 
dependence of the light induced width for the Je =1, 3, 5 
and 7 levels with the same vibrational quantum number as 
the E{Je = l)/h = -2235 MHz level. The light induced 
widths were calculated for scattering states with total angular 
momentum Jg = Jg — 1. The widths for the same scattering 
partial wave Jg but with Jg = Jg + 1 will be very similar in 
energy variation. The g-wave resonance can be clearly seen 
in upper panel of Fig. [3] 

To see how the photoassociation spectrum changes with the 
temperature of the gas sample we have calculated the ther- 
mally averaged light induced trap-loss coefficients for tem- 
peratures r =10 /iK, 100 /iK and 1 mK; see Fig. [3] These 
correspond to kT/h of 210 kHz, 2.1 MHz, and 21 MHz re- 
spectively. At temperature 10 fiK the spectrum is narrow, 
with a width determined mainly by the spontaneous emission 



rate, and is dominated by s-wave scattering, for which only 
the Je = 1 state is visible. The calculated spectrum at tem- 
perature T ~ 100 nK has three lines. One is the Je = 1 line 
due to s-wave scattering. This line has a normal thermal line 
shape 1I42I1 . Two other strong and sharp lines corresponding 
to Je = 3 and Je = 5 bound states are due to the g-wave 
resonance shown in upper pan el of Fig. [3] Both lines have 
subthermal width; see Ref. 157\. The shape of these two lines 
is determined mostly by the shape of the g-wave resonance. 
Finally, at a temperature of 1 mK one can clearly see the very 
broad thermally broadened lines coming from s and d-wave 
partial waves. The d-wave feature is much broader than the s- 
wave one. This is caused by the fact that for collision energy 
around 1 mK the light induced width for the s-wave decreases 
with an increase of the collision energy, while for the d-wave 
the light induced width increases. This affects the shape of 
the line and leads to some shifting and widening of the line 
coming from the d-wave compared to that coming from the 
s-wave. In addition, the g-wave shape resonance gives rise to 
a sharp subthermal feature. 

The second example is for ^^^Yb, which has a measured 
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FIG. 3: (Color online) Results for the excited 0^ bound states of the 
"^Yb2 molecule near the Je = 1 level at E/h = -2237 MHz; 
other rotational levels with the same vibrational quantum number 
are evident. The top panel shows contributions to the light induced 
width Veg coming from the optical coupling of the ground scatter- 



ing states having Jg 



la 



(s-wave), 2 (rf-wave), 4 (g-wave). 



and 6 (i-wave) with excited bound states having Je = Jg + 1, as a 
function of the collision energy £r, calculated for the light intensity 
1 W/cm'^. The s-wave scattering length is quite large in magnitude, 
—598(63) ao IT^ . so that this system has a g-wave shape resonance 
near er/fcs = 200 /iK that leads to resonantly enhanced J^ = 3 and 
5 features in the photoassociation spectrum. The lower three panels 
show contributions to the fully thermally averaged loss rate K com- 
ing from the optical coupling of the ground scattering states having 
Jg = ^9 ~ (s-wave), 2 (d-wave), 4 (p-wave), and 6 (i-wave) with 
excited bound states having Je = Jg + 1 and Je = Jg — 1, as a 
function of the laser detuning from the atomic resonance, calculated 
for various temperatures T = 10 ^K, 100 ^K, and 1 mK and the 
light intensity 0.1 mW/cm^. 



scattering length relatively close to a ITsl Issll . In such case 
one can expect that both a d-wave and an i-wave resonance 
can occur near threshold. As with Fig. [3] for ^^^Yb Fig. |4] 
shows the light induced width for the ^^"^ Yb bound-states with 
the total angular momentum Jg =1, 3, 5 and 7, with the same 
vibronic quantum number as the level at E{Je = l)/h = 
— 1973 MHz. The figure clearly shows the d and i-wave res- 
onances. As for ^^^Yb, the spectrum at a temperature of 10 
/iK is dominated by the s-wave component. However, unlike 
^^^Yb, the d-wave shape resonance leads to a weak Jg = 3 
feature even at such low temperature. The spectra at T = 100 
/iK and 1 mK clearly show the influence of the d-wave reso- 



FIG. 4: (Color online) Results for the excited OJ bound state of the 
^''''Yb2 molecule near the Je = 1 level at E/h = -1974 MHz. The 
panels are the same as described in the caption of Fig[3] The s-wave 
scattering length is 105(2) ao, which is about 1.4d, and there is a 
broad d-wave shape resonance near Er/ks ~ 220 /iK near the top 
of the d-wave centrifugal barrier IT^ . 



nance. This resonance had a crucial role in the correct inter- 
pretation of photoassociation spectra near the resonance tran- 
sition ^5*0 Pi and in the determination of the scattering 
length for this isotope Isst]. When the spectrum is dominated 
by the features due to the ground state d-wave, the intensity 
ratio of the Je = 3 and Je = 1 features should be 3/2 from Eq. 
(fl4] i. Such d-wave doublets in ^^"^Yb were observed by Tojo 
et al lfl3ll . who found that the transition to the bound state with 
Je = 3 can be stronger than the transition to Jg = 1 at rela- 
tively low collisions energies corresponding to a temperature 
of about 25 /iK, what can be explained by our model. 

The spectrum at temperature 1 mK is also dominated by 
the d-wave component. However one can notice sharp struc- 
ture on the top of a weak and broad feature due to the Jg = 5 
bound-state, mostly due to the ground state g-wave. The sharp 
structure is a consequence of the i-wave resonance. The i- 
wave component has a subthermal width connected to the i- 
wave resonance seen in upper panel of Fig. H] One can also 
notice a weak line corresponding to Je = 7 bound state sup- 
ported by this i-wave resonance. 

Figure|5]shows our calculations for Yb. This isotope has 
a relatively small negative scattering length ifTsll . The figure 
shows the light induced width for ^^^Yb bound-states with 
the total angular momentum Je =1, 3, 5 and 7, with the same 
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FIG. 5: (Color online) Results for the excited 0^ bound state of the 
"''Yb2 molecule near the Je = 1 level at E/h = -1739 MHz. 
The panels are the same as described in the caption of Fig [5] The 
s-wave scattering length is —24(5) ao, not near any special value, 
but there is a broad g-wave shape resonance at relatively high en- 
ergy near £r/fcs = 1 mK near the top of the p-wave centrifugal 
barrier 1131. 
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FIG. 6: (Color online) Results for the excited 0^ bound state of the 
"''Yb2 molecule near the Je = 1 level at E = -83 MHz. The 
panels are the same as described in the caption of Fig [3] In this case 
the binding energy is much smaller than in Figs. [3] ID and|5] and the 
Condon points are outside the centrifugal barriers except when the 
collision energy is very low; see Fig. [T] 



vibronic quantum number as the level at E{Je = l)/h = 
— 1737 MHz. Although the scattering length is not near any 
special value for van der Waals quantum defect theory, a wide 
g-wave resonance is clearly seen relatively high in energy near 
1 mK. The spectrum at 10 /iK is typical with only an s-wave 
line that can be observed. Upon increasing the temperature to 
100 /iK the spectrum is still dominated by the s-wave compo- 
nent, although noticeable d-wave components begin to appear 
Moreover very weak 5-wave components occur in this spec- 
trum. The picture changes dramatically when the temperature 
increases to 1 mK, where the g-wave components dominate 
the spectrum. 

Finally we would like to show how the spectrum is affected 
when the Condon point for the transition is moved outside the 
region of the centrifugal barrier. This is done by choosing a 
level with much smaller binding energy in our previous ex- 
amples. Figure |6] shows the light induced width for the ^^^Yb 
bound-states with the total angular momentum Jg =1, 3, 5 
and 7, with the same vibronic quantum number as the level at 
E{Je = l)/h = -83 MHz. The upper panel of Fig. |6]shows 
that there are no shape resonances, since there are no ampli- 
tude enhancements from being behind the centrifugal barrier 
The lower panels of Fig. |6]show the calculated spectra at 10 



fiK, 100 /iK, and 1 mK. The higher temperature 100 ^K spec- 
trum is dominated by broad d-wave components, on top of 
which a sharp subthermal s-wave component is clearly seen. 
This is because there is a node in the s-wave scattering wave 
function at the Condon point at relatively low scattering en- 
ergy, unlike for other partial waves. This sharp structure is 
a nice example of the subthermal line shapes discussed by 
Machholm et al. IstIi in the context of alkaline earth pho- 
toassociation. Increasing the temperature to 1 mK leads to a 
quasi-continuum as several partial waves contribute, as seen 
from Fig. |6] 

The examples in this section show that the photoassociation 
spectra of various isotopes of the same species can differ qual- 
itatively. This variety is directly connected with the quantum 
defect in the ground electronic state which is dependent on 
reduced mass and scattering length of the colliding species. 
This sensitivity of the rotational structure of photoassociation 
lines can help in the determination or verification of the scat- 
tering length, as it was recently done in case of calcium by 
Vogtef a/, il. 

A realistic description of near threshold interaction in the 
ground and excited electronic states of the Ybo molecule is 
possible thanks to the experimental data ifisi [l5ll and the cal- 
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FIG. 7: (Color online) The optical lengths l°g^ of the s-wave pho- 
toassociation features as a function of the detuning A of the reso- 
nance from the atomic transition for the several bosonic isotopes of 
Yb. The optical length determines the strength of an optical Fesh- 
bach resonance. The optical length varies linearly with laser inten- 
sity / and these values have been normalized to / = 1 W/cm^ for 

Er = 0. 



culations we have shown here. Our models of the interaction 
potentials can be used to predict the magnitude of the optical 
lengths in Eq. (|6|l for the set of optical Feshbach resonances 
coiTcsponding to the photoassociation lines of the different Yb 
bosonic isotopes (fermionic isotopes are more complicated 
because of hyperfine structure in the excited state 111411 ). Fig- 
ure |7] shows our calculated optical lengths for the homonu- 
clear bosonic isotopes of Yb. The isotope ^^^Yb offers the 
strongest optical Feshbach resonance because of its large neg- 
ative scattering length, which enhances the amplitude of the 
s-wave ground state wave function in the region of the Con- 
don points. Bound levels near the excited state threshold have 
optical lengths at 1 W cm~^ on the order of 10^ ap, similar to a 
value measured for a weakly bound excited level of **Sr IstIi . 
This suggests that these resonances may be of practical use for 
changing the scattering length for useful time scales while re- 
ducing spontaneous emission loss processes by using large de- 
tuning. Enomotoef a/. lfl6ll have experimentally demonstrated 
the possibility of some degree of optical Feshbach control in 
ultracold Yb gases. Our predicted optical lengths should be 
helpful in choosing good resonances for controlling ultracold 
Yb collisions by light. 



VIII. CONCLUSION 

One could naively think that this is a simple system so that 
not much difference should be observed between the spectra 
of three bosonic isotopes Hke ^ Yb, ^ ''^Yb, and ^''^Yb. How- 
ever, our calculations show that for temperatures about 100 
/LtK the spectra for various isotopes are qualitatively differ- 
ent. This difference is manifested by the fact that different 
rotational lines are apparent in different isotopes. Moreover 
the shapes of photoassociation lines are also affected. Some 



lines have subthermal widths depending on the isotope. These 
differences can be well understood as due to the isotopic vari- 
ation in the properties of the scattering wave function in the 
ground electronic state. There is a clear connection between 
shape resonances in ground scattering states and the appear- 
ance of strong lines coming from higher partial waves at rel- 
atively low temperature of the order 100 fiK. However, this 
only happens for excited bound states that have sufficiently 
short-range Condon points inside the location of the centrifu- 
gal baiTier. On the other hand, near-threshold bound states 
with long-range Condon points outside the centrifugal bar- 
rier do not show resonance enhancement. For example, we 
showed that the presence of Jg = 3 lines for ^^^Yb at small 
detunings does not indicate that a d-wave resonance occurs. 
Thus, a good experimental method to determine whether res- 
onances are present or not would be to measure photoasso- 
ciation hnes for more deeply bound excited states that have 
turning points near to or less that the ground state centrifu- 
gal barrier. The existence of resonances correlates with the 
approximate value of the s-wave scattering length. 



The intercombination photoassociation spectra of Yb are 
qualitatively different from spectra for group 11 atoms such as 
Ca and Sr. For the case of Ca theory shows that the excited 
state potential in the turning point range for levels near the 
dissociation limit is dominated by the van der Waals poten- 
tial, since the resonant dipolar interaction is so small ll22ll59ll . 
The reflection approximation is not applicable in such a case 
lliUl . Strontium l37ll is a system which is half way between 
Ca and Yb. Only the last three Sr bound states closest to the 
dissociation limit can be treated as dominated by the resonant 
dipole interaction, while more deeply bound levels are deter- 
mined mostly by the van der Waals part of the interaction. 
Photoassociation can be observed in Ca or Sr for two series of 
bound states of 0+ and 1„ symmetry, respectively. Although 
the resonant dipole interaction for the lu state is repulsive, its 
weakness allows it to be overcome by the attractive van der 
Waals interaction. In contrast the Yb interaction in the excited 
electronic state with dissociation limit ^Sa Pi is dominated 
by the resonance interaction, so that only one potential curve 
of 0+ symmetry is attractive at long range and supports a se- 
ries of detectable bound states. The reflection approximation 
is well applicable in the case of 0+ levels of Yb. Since the 1„ 
potential becomes attractive at much shorter range of the inter- 
atomic separation, the density of bound states near threshold 
is much less, and only one bound state with 1„ symmetry is 
expected near within a few GHz of threshold. 



The heavier group lib elements such as Cd and Hg will be 
similar to Yb. For both species the interaction in the excited 
state is dominated by the resonance interaction, and therefore 
one can expect that the reflection approximation will be appli- 
cable. Moreover both species have numerous stable isotopes. 
Consequently, mass tuning of the scattering length should be 
applicable for both of these species, and the results obtained 
in this work can be treated as universal and generic for Cd and 
Hg. 
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